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Abstract: We present the automation of a subtraction method for photon radiation 
using the dipole formalism within the MadGraph framework. The subtraction terms are 
implemented both in dimensional regularization and mass regularization for massless and 
massive cases and non-collinear-safe observables are accounted for. 
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1. Introduction 

A precise theoretical description of cross sections at particle colliders is required for the 
reliable extraction of fundamental parameters, for indirect searches for new physics effects, 
and for a well-founded assessment of background reactions to new physics signatures and 
their uncertainties. This precision is achieved by inclusion of higher order perturbative 
corrections from QCD and from the electroweak theory. At hadron colliders, especially the 
inclusion of next-to-leading order (NLO) QCD appears to be mandatory for many studies. 
Much effort has been made in the recent past to automate NLO QCD calculations (see [1,2] 
for reviews). 

A priori, NLO electroweak corrections can be expected to be smaller than QCD cor- 
rections, owing to the smaller coupling constant. In a variety of situations, this naive argu- 
mentation does however not apply. In QCD, quarks and gluons can not be distinguished, 
and collinear gluon radiation off an outgoing parton can not be identified experimentally. 
Collinear photon radiation off a lepton leaves a signature in the detector and can thus be 
identified. As a result, the cancellation of photonic contributions in the virtual and real 
collinear contributions at NLO is incomplete, and the observable is non-collinear-safe [3]. 
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The final state photon selection cuts must be implemented consistently in the theoretical 
description, both at fixed order and in a parton shower framework [4]. Depending on the 
treatment of photon radiation in the experimental analysis, the magnitude and shape of 
the electroweak corrections can differ substantially, due to the occurrence of large loga- 
rithmic terms in the photon cuts, as discussed in detail for instance in the electroweak 
corrections to gauge boson production and decay [5,6]. At centre-of-mass energies above 
the weak gauge boson masses, similar incomplete cancellations occur in the genuine weak 
corrections, giving rise to large logarithms in the gauge boson masses. 

NLO electroweak corrections for hadron collider observables have up to now been de- 
rived for gauge boson production [5,6], for gauge boson pair production [7], for Higgs 
boson production and decay [8-13], for top quark pair production [14, 15], single top pro- 
duction [16] and for supersymmetric particle production [17-19]. Much effort has been 
made recently towards electroweak corrections to final states including jets. For these 
final states, the experimental treatment of photons observed inside jets [20] can yield non- 
collinear-safe configurations, resulting in potentially large corrections. A first calculation of 
this type was accomplished recently with the NLO corrections to T4 /± +jet final states [21]. 

NLO electroweak corrections require two ingredients: the one-loop virtual corrections 
and the real radiation corrections with one extra final state particle. While the one-loop vir- 
tual corrections contain explicit infrared poles from the loop integration, the real radiation 
contributions develop infrared singularities only when integrated over the real radiation 
phase space. To handle these real radiation singularities in a process-independent manner, 
subtraction formalisms have been developed initially in the context of QCD calculations: 
residue (or FKS) subtraction [22] and variants thereof [23], dipole subtraction [24,25] and 
antenna subtraction [26-29], all within dimensional regularization. These can in principle 
be carried over to electroweak calculations. However, some care has to be taken since most 
electroweak calculations use mass regularization. The dipole subtraction method has been 
extended to include photonic radiation in [30]. 

Since NLO electroweak corrections are required for a substantial number of processes, 
their automation is very much desirable. Driven mainly by QCD calculations, the au- 
tomation of one-loop corrections for one-loop multi-parton amplitudes has seen enormous 
progress in the recent past, and first fully automated programs for the calculation of one- 
loop multi-parton amplitudes are becoming available with the packages CutTools [31,32], 
BlackHat [33], Rocket [34], GOLEM [35] and Samurai [36], as well as independent li- 
braries [37]. 

Automation of the real radiation contributions is normally based on a leading order 
parton- level event generator framework such as MadGraph/MadEvent [38-40], ALPGEN 
[41], CompHEP/CalcHEP [42]/ [43], SHERPA [44,45], HELAC [46] or WHIZARD [47]. 
The generation of dipole terms for subtracting the singular behaviour from the real ra- 
diation subprocesses has been automated in various event generators: in the SHERPA 
framework [48], the TeVJet framework [49], the HELAC framework [50] and in the form of 
independent libraries [51,52] interfaced to MadGraph. The MadDipole package [53] pro- 
vides an implementation within MadGraph. An implementation of the residue subtraction 
method is also available within MadGraph [54]. For a full NLO calculation, the dipole 
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terms have to be integrated over the dipole phase space, and added with the virtual correc- 
tions to obtain the cancellation of infrared singularities. The HELAC implementation [50] 
and the MadDipole package [55] provide these integrated dipole terms including all masses 
and possible phase-space restrictions. 

It is the purpose of this work to extend the MadDipole package towards electroweak 
NLO corrections by automatically generating the unintegrated and integrated dipole sub- 
traction terms for photonic radiation, as defined in [30], and including non-collinear-safe 
configurations [3]. 

This paper is structured as follows: in Section ^J, we review the differences between mass 
regularization and dimensional regularization, and their implications for the dipole terms. 
The treatment of photonic effects in the parton distributions is described in Section ||[ 
and the dipole terms for non-collinear-safe observables are introduced in Section ||. The 
implementation of photon radiation in MadDipole is documented in Section [|, and checks 
of it are described in Section ||. 

2. Dimensional regularization vs. mass regularization 

Infrared singularities are an inherent property of quantum field theories. In the dipole 
formalism this becomes exlicit when integrating over the one particle phase space to obtain 
the integrated subtraction terms. As the integration over this phase space is divergent one 
needs a prescription to parametrize the singularities. For QCD calculations dimensional 
regularization is the most widely used method. The singularities are present in poles of 
the form 1/e 2 and 1/e. Electroweak calculations are mostly done using mass regularization 
where the actual masses of the fermions and a hypothetical mass parameter of the photon 
are used as infrared cut-offs. In this case the singularities show up as logarithms of the 
mass parameter. 

In our implementation we offer both possibilities and the user can select which method is 
used. The implementation within dimensional regularization works in the same way as for 
the QCD version of MadDipole, described in [53,55]. 

In handling photon radiation, the main differences are in the initial state radiation collinear 
singularities, which lack counterparts in the virtual corrections. In QCD these singularities 
are absorbed by a redefinition of the parton distribution functions. Integrated subtraction 
terms for processes with initial state partons contain distributions in the momentum frac- 
tion x of the parton entering the hard scattering process (for details see [24,25]). A priori 
the reduced matrix element also depends on x, a feature that bears disadvantages for the 
numerical implementation: the matrix element has to be evaluated at x ^ 1 and at x = 1. 
In the presence of a parton distribution function this problem can be circumvented by a 
transformation of variables, thus shifting the dependence of this variable into the parton 
distribution function. The details and how this is done in MadDipole is described in [55]. 

The major new aspects in the implementation of photon radiation are the following: 

1. Initial state radiation singularities in electroweak processes are regulated by the masses 

of leptons entering the hard process, resulting in potentially large mass logarithms. In 
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contrast to the divergencies absorbed into the parton distribution functions, these mass 
logarithms have physical significance. Consequently, for electroweak corrections to pro- 
cesses with leptons in the initial state, the user has to choose mass regularization as the 
appropriate scheme. Furthermore, it is not possible to remove the x-dependence from the 
hard scattering matrix element by a variable redefinition, which thus has to be evaluated 
up to three times: at x = 1, and at the values of x corresponding to the initial state 
radiation off each of the two incoming momenta. 

2. For configurations containing only a collinear but not a soft divergence, the summation 
over the dipole spectator momenta is obsolete, since these configurations can be accounted 
for by a single dipole term. In the QCD MadDipole implementation, this summation is 
nevertheless carried out in order to facilitate the colour management. In the electroweak 
MadDipole implementation, this summation is omitted. For the convenience of the user, 
all possible spectator terms are generated, and can be invoked manually. 



3. QED corrections to parton distribution functions 

For processes with incoming partons, initial state collinear singularities are absorbed into 
the parton distribution functions (PDFs). In dimensional regularization this collinear sin- 
gularity is of the form 

(4vr/i 2 ) e 1 



r(l - e) e 



Pji(x) , (3.1) 



where Pji denotes the Altarelli-Parisi splitting function for the j —> i splitting. 

For example, if one considers gluon radiation off an incoming quark, the absorption of the 



collinear singularity of |3.1| into the quark PDF leads to 
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(3.2) 



where hf is the factorization scale. The value of the coefficient C qq ( ~J defines the sub- 
traction scheme. The MS scheme which is mainly used in QCD just absorbs the singular 
terms and some trivial constants, so 

Cjp(x) = . (3.3) 

For the subtraction of photon radiation, the equivalent expressions for QED corrections 
using mass regularization are needed. This can be derived by calculating the integrated 
subtraction terms both in dimensional and mass regularization and comparing the results. 
A detailed derivation can be found in [8,56,57], so we only quote the results here. According 
to Eqns. (3.19) and (3.20) of [56], to include inital state radiation at NLO QED of a parton a 
out of a hadron h with longitudinal momentum fraction x, the parton distribution functions 
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are modified by 

x {ln(4) lP„(z)} + - lP„(z) (21n(l - z) + 1)] + + Cjf (*) } 

a=q,q x 

X { ln 0f) P 7/W-^7/W(21n^+l) + C 7 D / IS (z) } , 

with the splitting functions 

P ff^ = TT7> P fl {z)=z 2 + (l-zf, P lf {z)= l + {l ~ Z? . (3.6) 

The only set of PDFs which include both QCD and QED NLO corrections is the MRST2004QED 
[58]. As has been pointed out in [57], a consistent use of this set of PDFs requires the use 
of the DIS scheme for QED calculations. The coefficient functions for the DIS scheme are 
given by 



, , ,'l-z\ 3\ 9 + 5z 



c ff i*) : z j 4 y • 4 

C?™(z) = P h (z) ln (^) - ^z 2 + 8z - 1 , 

C»?{z) = -C?f(z) . (3.7) 

It should be noted that the integrated subtraction terms as they appear in MadDipole are 
different for processes with incoming PDF and for processes without. The introduction of 
the PDF does not only affect the logarithmically divergent terms, where the mass of the 
particle is replaced by the factorization scale, but also finite parts are different due to the 
coefficient functions present in the PDFs. Whether a PDF is present or not is determined 
automatically by MadDipole in the case of leptons or QCD particles. Any QCD particle 
(quark, anti-quark, gluon) is assumed to be a constituent of a hadron whereas this is 
certainly not true for leptons. A special case is the incoming photon, because it may come 
out of a hadron but could as well be a free particle. Here the user has to choose whether 
to use a PDF or not (see section |5|) . 

4. Non-collinear-safe observables 

Final state collinear photons can often be distinguished from other final state particles. 
For example, the energy fractions in a collinear muon-photon system can be disentangled 
by comparing calorimetric and tracking information, and highly energetic photons inside 
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QCD jets can be identified and vetoed. Therefore, in electroweak calculations a collinear 
photon is not necessarily treated fully inclusively [3]. Consequently, the integration range 
within the cone of a photon radiated off a charged particle is not the full range of the 
momentum fraction of the charged particle. This momentum fraction is parametrized by 
the variable Zij for a final-final constellation and Zi a for a final-inital one. Treating the 
photon not fully inclusively is now equivalent to imposing a cut on this variable. In this 
case, the collinear singularity associated with the photon emission off a massless fermion is 
not cancelled, rendering the observable non-collinear-safe. In analogy to the factorization 
of initial state collinear singularities into parton distributions, this singularity is absorbed 
into a photon fragmentation function [20,59]. An alternative photon isolation procedure 
has been proposed in the form of a dynamical cone-size [60], which would in principle shield 
the collinear singularity and eliminate the need for the photon fragmentation function. In 
practise, this criterion turns out to be difficult to implement for a finite-resolution detector. 
The treatment of fragmentation processes in the dipole formalism was first developed in [24] 
for dimensional regularization, and extended to mass regularization in [3]. We largely fol- 
low the latter paper for notation and conventions and refer to this reference for further 
details. Here we collect those results which were important for our implementation. 

For observables with identified final state photons, the real emission matrix element and 
the unintegrated subtraction terms are multiplied by a step function that cuts off the 
integration at a given value for z. One therefore has 



The integrated dipole terms are thus not integrated over z, but remain functions of this 
variable. 

The MadDipole implementation usually allows for the introduction of a phase space cut, a, 
to ensure that dipoles are subtracted only in the vicinity of the singular regions. In our im- 
plementation, non-collinear-safe observables can only be calculated with all a-parameters 
set to one, since including the ©-function related to the a-parameter in them would lead 
to very involved phase space constraints, which prevent the analytical integration of the 
subtraction terms. 

For non-collinear safe-observables, the emitter is always in the final state, while the 
spectator can be in final or initial state. Both cases were derived in [3] and are summarized 
below. 
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4.1 final-final 

In the final-final case the integrated splitting function in mass regularization is given by 



F 4 



-I 



V2(zij) 



dyij (l - yij) g^ U r\pi,Pj, k) . 



(4.2) 



For an appropriate numerical treatment the soft singularity at z = 1 is regularized by 
introducing a plus distribution for z, 



-y(sub) 
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(4.3) 



The integration boundaries can be derived from the phase space, for details see Appendix 
B of [30] . In the limit mj — > and rrij , m 7 = one finds 

mf (1 - z) 



Vi{z) = 

Performing the integration leads to 

g^i\p^z) = p n {z) 



P 2 z 

IJ 



V2(Z) = 1. 



(4.4) 



In 



4f 



- 1 



+ (l + z)ln(l-z) 



with 



1 + z 2 
1-z 



(4.5) 



(4.6) 



The case where rrij is not zero is more involved however there is no conceptual difference 
to the case of a massless spectator so we do not quote the result here. It can be found in 
Appendix A of [3]. 

The derivation of the corresponding expressions in dimensional regularization starts 
from the dipole splitting function 



< V ijjk (zi;y) > 



which integrates to 



(1 - Vij,k) L 1 - z i( l ~ Vij,k) 



- e P ff (z) + P ff (z) Hz) + ln(l - z) + (1 - z) 



V(e,z) 

In the case where the spectator is massive, the splitting function takes the form 



< Vij,k(zj]y) > 
integrating to 



1 



2 ^ij,k 

Rij(y) [1-Zj(l -y ijjk ) v ij:k 



(1 + Zj + e(l - Zj)) 



(4.7) 



(4.t 



(4.9) 



- {z 2 (- log (8(ci - c 2 )) + log(4(c(c + c 2 ) - ciy 2 )) + log(c 3 )) 

+ log (ci - c 2 ) - log(c(c + c 2 ) - ciy 2 ) - 2 (z 2 - l) log(c) - log(c 3 ) 
+21og((y 2 - l)z + 1) - 21og(l - z) + log(2)} + P ff (z) log(-cy 2 (z - l)z) -z + 1, 

(4.10) 
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where we used the abbreviations 



c=1 -mL c i = 1 -/"fc> c 2 = \Jc 2 (yl + 1) - 2ciy 2 , c 3 = ci - c(c y 2 + c 2 ) (4.11) 
and 



to compactify the result. 
4.2 final-initial 



4 



2z(l - z) 



(4.12) 



For an initial state spectator the subtraction terms depend on the variable X{ a and are 
divergent for Xi a — > 1. This is taken into account by splitting off the divergent parts using 
the (+)-distribution. The correction terms for non-collinear-safe observables now come 
along with an additional distribution in z so the the full result reads as a combination of 
the two. 



p2 fX\ 

ia 

2 



f Xl r z ^ x ) ( snh \ 

dx dzg™ T (pi,p a ,k) 

J0 Jz 1 (x) 

= f dx f • 
Jo Jo 



dziG^\Pl)6(l-x)S(l-z)+ glZ°>(P? a ,x) 5(1 -z) 



;(sub)/ p 2 



with 



_(sub) / x 
9ia,+ ( x > z ) 



+ 



1 



_(sub) 



(x,z) 



(x,z) 



(4.13) 
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+ (l + z)ln(l-«), 



G^(Pl,z) = l-z. 



(4.14) 



The first contribution, gf^\x,z) is identical in mass regularization and dimensional 
regularization. 

To evaluate the equivalent expression of Qfar\Pfa^ z ) in dimensional regularization we use 
the splitting function and phase space given in Eqns. (5.39) and (5.48) of [24] and find 



PffW 



z 2 + 3 



,log(2-z) 



V(e, z) = l^-L + P ff (z) log(z) + - T — log(l -z)- 2 ^— ' + 1 

4.3 Factorization of the collinear final state singularity 



— z . 



(4.15) 



The integrated dipole factors contain a collinear singularity. For fully inclusive observables, 
the dipole factors would be integrated over the photon momentum fraction z and added 
with the virtual photon corrections, rendering the result finite. Since non-collinear-safe 
observables explicitly require (or veto) the final state photon, these cancellations no longer 
take place. In the case of photon radiation off massive fermions, this results in a large 
logarithm in the fermion mass, which is present in the above mass-regularization results. 
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For massless fermions, especially for light quarks, this turns into a collinear singularity, 
and one has to include a contribution from the fragmentation of a final state fermion into a 
photon-fermion system of low invariant mass. The quark-to-photon fragmentation function 
(which depends on the photon momentum fraction (1 — z)) is an empirically determined 
quantity [20,59,61,62]. It contains a counterterm for the collinear singularity, whose explicit 
form depends on the regularization scheme. The counterterm is accounted for by adding 
the bare fragmentation function to the integrated dipole terms. In mass regularization, the 
bare fragmentation function reads [63]: 



D 



bare, MR/ 
<3->7 * 



1 - Z ) = D q -^. 7 (1 - Z, flp) + 



aQl 




+ 1 



(4.16) 



while the corresponding expression in dimensional regularization is [20]: 

">2 -, / a-.,2\ e 



D 



bare,DR 
9^7 



(1 - z) = D q ^ 7 (l - z, fj, F ) + 



1 



aQ z q 1 /4vr^ 2 
~2tT1 V ti F J r(i -e 



-P ff (z). (4.17) 



It can be seen that the sum of the integrated dipole terms and the bare fragmentation func- 
tion is independent on the regularization scheme. The scale-dependent non-perturbative 
fragmentation function D q ^(l — z,fip) has to be determined from data. 
Our implementation offers the possibility to include the fragmentation function. For details 
see section || 



5. Implementation and how to use it 

The installation and running of the new package is identical to the QCD version: 

1. Download the MadDipole package (version 4.4.35 or later), MG_ME_DIP_V4 . 4 . ?? . tar . gz, 
from one of the MadGraph websites, e.g., http://madgraph.hep.uiuc.edu/. 

2. Extract and run make in the MadGraphll directory. 

3. Copy the Template directory into a new directory, e.g., MyProcDir to ensure that 
you always have a clean copy of the Template directory. 

4. Go to the new MyProcDir directory and specify your process in the file . /Cards/proc_card . dat. 
This is the (n + l)-particle process you require the subtraction term for. In this file 

there is also a variabe DipolePhotonPdf which is set to TRUE by default. If so it is 
assumed that initial state photons are coming from a hadron and therefore a PDF is 
used. 

5. Running . /bin/newprocess_qed generates the code for the (n + l)-particle matrix 
element and for all dipole terms and their integrated versions. After running this 
you will find a newly generated directory . /SubProcesses/PCLyourprocess (e.g., 
. /SubProcesses/PO_e+e-_uuxa) which contains all required files. 
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In the . /SubProcesses/PCLyourprocess directory all the files relevant to that par- 
ticular subprocesses are generated. In particular this includes the (n + 1) particle matrix 
elements in the file matrix . f and the dipoles in the files dipolqed??? . f, where ??? stands 
for a number starting from 001. Furthermore the directory has two files, dipolsumqed. f 
and intdipolesqed.f , where the sum of the dipoles and their integrated versions are cal- 
culated, respectively. 

The routines contained in intdipolesqed.f which calculate the integrated subtraction 
terms have a slightly different syntax compared to the QCD version so we give some more 
details here. As for QCD the file contains two routines 

intdipolesqed(P, PI , P2 , X, Z, Zi, PSWGT EPSSQ, EPS, FIN) and 
intdipolesqedfinite(P, PI , P2 , X, Z, PSWGT EPSSQ , EPS, FIN). 

The difference is that the routines need three phase space points as input variables. This is 
due to the problem explained in section 1 , that in the absence of a PDF the reduced matrix 
element depends on the longitudinal momentum fraction x (here named as Z). Therefore 
the phase space point P(0 : 3,nexternal) denotes the point in an n-particle phase space 
with x = 1. The next variable, Pl(0:3,nexternal) denotes the phase space point with an 
initial state radiation from the first incoming particle, which means its momentum is given 
by x ■ pi where p\ is the full first momentum as it appears in the first phase space point. 
The third phase space point is the same as second but with the radiation off the second 
incoming particle. 

These phase space points have to be provided by the user. In the case where there are 
PDFs present for both incoming particles we use the same method as in QCD, namely 
shifting the dependence of the matrix element on the momentum fraction x into the PDF. 
In that case only the first phase space point P(0:3,nexternal) has to be provided, the 
others are not used. 

The second difference to QCD is the appearence of the variable Zi. It is only needed 
for non-collinear-safe observables (see section ||) and denotes the integration variable z as 
defined in that context. We assume that the integration over Zi is carried out numerically 
and therefore this variable has to be provided by the user. 

Several important parameters are included in the file dipole . inc which need to be set 
by the user. 

The first parameter decides which regularization scheme should be used. It can be either 
DREG or MREG which stands for dimensional regularization and mass regularization respec- 
tively. Note that in the case where the reduced matrix element contains leptons in the 
initial state, one has to choose mass regularization. 

The second parameter is a logical variable photonpdf which determines whether a PDF 
should be used for initial state photons. Whether this parameter is true or false per default 
is determined by the choice the user made in param_card.dat. If the according variable 
in the parameter card is set to true, this variable in dipole . inc will be set to true au- 
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tomatically when the code is generated. So for consistency this parameter should not be 
changed by the user. The next variable is a real array massint which contains the particle 
masses and which have to be given by the user. More precisely, what is required are the 
masses of the particles as they appear in the integrated subtraction terms. As an example 
we can take e + e~ — > tij. For the real emission matrix element and for the unintegrated 
dipoles one would usually treat the electron as a massless particle whereas the top-quark 
would be treated as a massive particle. The integrated subtraction terms however require 
an electron mass to regulate the infrared singularities. These regulator masses are declared 
in the array. So in this example one would fill the first two entries with the electron mass, 
and the next two entries would be the top mass. A photon should be treated massless in 
this array, so the last entry would be zero. All units are in GeV. Per default all entries are 
set to one to ensure a finite result. 

The next parameter is the photon mass parameter, which is also needed to regulate log- 
arithmic divergencies. By default it is set to unity so that all logarithms containing the 
photon mass vanish. 

The following parameter is a logical variable ncs. If set to true it calculates the correction 
terms for non-collinear safe observables. This requires of course that the user has to provide 
appropriate values for the variable Zi described above. 

The last parameter is given by a logical variable phot onf rag which is set to true by default. 
If set to true the second terms on the right hand side of Eqns. ( 4.16|/4~T^ ) are added to 
the integrated dipoles and a template function is called which can be used by the user to 
provide the non-perturbative parts of the fragmentation function, given by the first term 
on the right hand side of Eqns. ( f4.16 / 4~T^ ), This template function is described in section 

o. 



5.1 Implementation of non-collinear-safe case 

The calculation of non-collinear-safe observable requires to introduce additional cuts which 



act on the matrix element as well as on the subtraction terms, as in Eq. (4.1). In MadDipole 
these cuts can be imposed using the predefined cut routines, no additional routines are 
introduced. The correct handling of the G-function has to be done by the user. 
Also the integrated subtraction terms receive additional cuts in form of a 0-function. For 
the final-final case, the integrated subtraction terms are therefore of the form 



J d$i |A4 ub ,ij($i;>s)| 2 = -^Qi<JiQjOj J d$ ,ij J 



dz 



z) + 



'tJ,T 



(5.1) 



X \M (pi,Pj;TKi)\ 2 @ cut (pi = Zpi,k = (1 - z)Pi,Pj,{krSj: 



The parts that are z-dependent are the (+)-distribution and the 0-function. This means 
that this 0-function regularizes the singularity contained in the (+)-distribution at z = 1. 
We therefore need to evaluate the step function at two points, z ^ 1 and 2 = 1. As this can 
not be done within the predefined cut routines, we introduce an additional logical variable 
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• thetacut (P,zi) 



which can be either true or false. Its arguments are the phase space point and the variable 
z. As MadDipole uses a variable z already to denote the momentum fraction in the context 
of initial state radiation we rename this variable to zi. In the same way as for the 'normal' 
cut routine MadDipole provides an empty dummy routine in which the user can specify 
the desired cuts. MadDipole then takes care that it is evaluated at the two different points 
and combined with the subtraction terms according to (|5.1| ). The template of this function 
is included in the file intdipolesqed.f . 

In the same file we also provide a template that can be used for calculating the non- 
perturbative parts of the quark-to-photon fragmentation function as described in the pre- 
vious section. The template is given by the real function 

• DFRAG(zi) 

where the only input variable is z. Here the user should provide a call to the set of 
fragmentation functions of his choice. The result will then be added automatically to the 
integrated dipole terms. Per default, this routine returns zero so no non-perturbative terms 
are added. 



6. Checks 



Several checks have been performed to test the correctness of the different parts of the 
implementation. A first check of the unintegrated dipoles is of course to check the limits, 
i.e. the fact that the ratio of real emission matrix element over the sum of all subtraction 
terms, 

1 (6.1) 

when one approaches a soft- and/or collinear limit. This could easily be done using the 
already existing check program which is part of MadDipole (see [53] for details). 
This alone is not a sufficient test. Furthermore we compared results for single phase space 
points against existing results for the processes e + e~ —> 4/ with RacoonWW [64-66], 
e~7 — > e~fi + fi~ [3] and e + e~ — > tiH [67,68]. Thus all different splittings could be tested. 
This has been done for both the unintegrated dipoles as well as for the integrated ones. 
The aim of this package is to provide the necessary subtraction terms for single phase 
space points without the phase space integration. In particular the user is not restricted 
to MadEvent for the integration but can use MadDipole as a stand alone package and 
implement it in a completely different framework. A pointwise agreement with existing 
results should therefore be considered sufficient to test the implementation. The calculation 
of complete cross sections (with virtual corrections from other packages) is beyond the scope 
of this work. 

Furthermore we repeated the checks on the dependence of the cut parameter a as described 
in [55]. Although already being checked in the QCD version of MadDipole additional checks 
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(a) (b) 
Figure 1: Dependence of the cross section on the cut parameter a. All four a-parameters (final- 
final, final-initial, initial-final, initial-initial) are varied simultaneously. The dashed line denotes 
the real emission matrix element where the dipoles have been subtracted. The dotted line is the 
contribution from the finite parts of the integrated dipoles. The solid red line is the sum of the two 
contributions. As this is by no means a physical meaningful quantity we do not give actual numbers 
for the cross section. For the sum, the error bars of the Monte Carlo integration are added. 

were necessary because of two reasons. First of all the final- final dipole splitting function 
/* — > /7 is conventionally normalized differently than in QCD, which made it necessary 
to calculate the correction terms (see appendix). On the other hand having leptons in 
the initial state made it necessary to incorporate the concept of calculating amplitudes at 
different points [x = 1 and x ^ 1) which is not needed for QCD processes. Checking the 
independence of the total result on the phase space cut parameter provides a valuable test 
for this part of the implementation. 

In Fig.|| we show the dependence of the a-parameter for two processes, e + e~ — > uuj and 
e + e~ — > tty. These two processes contain the final-final splittings different from QCD and 
they also contain initial state radiation off leptons. What the plots show is the cross section 
as a function of the cut parameter. Actually, each splitting (final-final, final-initial, initial- 
final and initial-initial) has its own a-parameter which can be adjusted independently from 
each other. In this plot we change them simultaneously so that their values are always 
equal. 

The dashed line shows the cross section of the real emission matrix element where the 
dipoles have been subtracted so that the difference is finite. The dotted line shows the 
cross section of the finite parts of the integrated dipoles. Summing both contributions 
leads to the solid red line which shows an independence on the cut parameter to a very 
good accuracy over several orders of magnitude. 

7. Conclusions 

The MadDipole package automatically generates unintegrated and integrated subtraction 
terms required for NLO calculations. Previously, it provided all dipole subtraction terms 
for QCD radiation [53,55]. In this paper, we extended MadDipole to include photonic 
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radiation. The implementation is provided both in dimensional regularization and mass 
regularization. We include counterterms for the QED corrections to the parton distribu- 
tion functions, and allow for subtraction terms for non-collinear-safe observables, i.e. with 
identified final state photons. 

With the photonic radiation now included in MadDipole, the MadGraph/MadEvent 
framework can now be used for electroweak real radiation corrections to collider observables. 
This development provides an important step towards the full automation of electroweak 
NLO corrections. 
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A. Final-final /—>■/ + 7 dipole in dimensional regularization 

The final-final dipole factor for the radiation of a photon off a fermion is conventionally 
normalized differently from the corresponding QCD case. It is therefore documented in 
detail here. For the case where a photon is radiated from a massless fermion the auxiliary 
function is given by (Eqn. (3.1) of [30]) 



1 - Zij(l - yij) iJ y 

In this notation k is the momentum of the photon, pi and pj are the momenta of emitter 
and spectator. The variables and correspond to the variables y and z of the QCD 
case. Note that there is an additional factor of j^j- present in this case. 
The integrated term in dimensional regularization is obtained along the lines of the Catani- 
Seymour approach [24]. 

We want to use dimensional regluarisation therefore we are now following the Catani- 
Seymour approach. In that language the integrated term is given by Eqns. (5.27,5.28) 
of [24]: 

V,j, = I [*,«,,»)] ^ < V ti , t > s {^C)V,M ■ (A.2) 

and 

v ij (e) = J o dziMi-Zi)) y o -(i-y) V 87rag ^ ■ (A-3) 
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The transition can be made by replacing 

< V ijjk (zi;y) > _^ 1 



1 - Zij - e(l - Zij) 



(A.4) 



87ra S M 2e (1 - yij) |_1 - Zij{\ - yij) 

After integration and expanding in e we find 

^)=? + I + G-t) +o ( £1 )- < a - 5 > 

The massive case is more involved. We decompose the integral into an eikonal and a 
collinear part, 

Ijf,k{^j,^k;e) = 2I eik (/Xj,/x fe ;e) + I^ k (fj,j , fj, k ; e) , (A. 6) 

where we define the eikonal integral in analogy to Eq. (5.26) of [25] 

The difference to the CS dipoles is given by the additional factor Rij(y), which is defined 
in Eq. (4.3) of [30] 



Rij(y) = * —7=^ — , (A.8) 



with 

Pij =Pi+ Pj + k, Pfj = Pfj - mf -rrij- m 2 v \j = A(P? , mf , mf) . (A. 9) 
Performing the integration we find for the eikonal integral 

T eik (n ■ A — 1 lnp i xl/li ( {cj-v+\){c k -v+\) \ j ■ ( (c j +c k +2)(c k -y + + l)-C2 \ 

1 W],flk,e)- e $ "I" y |^12 \{ Cj +c k +2){c k -y + + l)-C2 ) + ^ \ (c+v+l)(c k +v+l) J 

+2 ln(p) ln((c, + c k + 2){c k - y + + 1) - c 2 ) + \ In* ( ^-S^^-l^ ) 

+ In(2S) ^ ( (cj+Cfc+ ffi-r i) - Ca ) + (ln( Cj - v + 1) In (a±p) 
- In (25( Cj -5 + 1)) ln(c fc - t; + 1) - In (2v) ln( Cj - 5 + 1)) + Li 2 (^gq^) 

+ Li 2 (^w) + + 5 + !) ln (f ) + 21n (^') hfe + y+ + 1) 
+^ " ^ (^Sr) " ^ ~ 2 m fo) ln(<* - y + + 1) 

+ 1 Mp J )ln( g||S^ ) -^n 2 (^g) -Li^ + Li,^) 

-Li 2 (pf ) - £ + Li 2 (-^TI)) - Info) In } + O(c) , (A.10) 

where we used the following abbreviations: 

= y/xil^lnl) ,p n ( Mj ., Mfc ) = yEg±^|Ep^ (n = j,/c), 



P = A/y^f, (A.H) 
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and 

y + = i - ; 2 "f-? k) 2 (A.12) 

is the upper limit of the y integration. Furthermore we define for brevity 



,2 _ ,,2 ' Cfc i ,,2 ,,2 ' 



c fc = 1 c 2 = J k . (A.13) 



The contribution from the collinear integral is simpler and we find 

I$ k (fi f ,n k ; e) = \ + \ (-4 log + 4~ 2 ^ + l) 



^ + 3^ + log(^(l - /x fe )) + O(e) . (A.14) 



(Mfc 

The results in dimensional regularization can be compared with the results in mass regu- 
larization given in Eq. (4.10) of [30]. Comparing the two results numerically we find full 
agreement. 



As the limit e — > does not commute with the limit of vanishing masses we also give 
the explicit result for the case when one mass is set to zero. In the case of a vanishing 
emitter mass we find 

/ cik (0, ^ e) = & - + 2Li 2 (jjjJLj) + In - 1) m(l - „ k ) 

+ ln(/i fc + l)(41n(l - Mfc ) - 31n( R .) + 21n( R , + 1)) - 7 -g + 0(e) , (A.15) 
whereas for vanishing spectator mass we find 

/ eik (^-, 0; e) = ^ + 2Li 2 ^1 - - Li 2 (^ 2 ) - 41og (l - ^) log(^) 

+ 31og 2 (/i,) + ^ + 0(e) . (A.16) 
The corresponding terms from the collinear integral read as 

l*f k (0,fi k ;e) = A + 1 - 31og(l - W ) + 0(e) , (A.17) 

and 

I*f k ( H , 0; e) = i + \ {ii) - 4 log (l - M 2 ) + 21og(/x j ) + 3} + 0(e) . (A.18) 
A.l a-cut parameter 

The idea and the use of the a-parameter is described in detail in [69,70]. As the final-final 
splitting function for the splitting /* — > fj is different from the one in QCD the results 
were not present in MadDipole so far. 

These correction terms only affect finite terms and have no influence on singular terms 
which means they can be calculated either in the mass regularization approach or in the 
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dimensional regularization approach. The correction terms are the same so one can use 
the method that appears to be more suitable. 

The general form of the correction terms in the language of dimensional regularization is 
given by 

I 7fAVj^k]t,a) = Iyf :k (Vj,Vk;e) + AI~ / f,k(L L jif i k;a) . (A. 19) 

The massless case is straightforward and we find 

A/ 7/)fc (0, 0; a) = -2Li 2 (l - a) - log 2 (a) - . (A.20) 

In the massive case we calculate eikonal integral and collinear integral seperately. For the 
case where both emitter and spectator are massive we find 



reik/ \ 1 for- ( v +c + cc k + c k - (c + c k + 2)y a + l 

2v [ \ (c + v + l)(cfc + v + 1) 

_ 2T . ( v 2 + c + cc k + c k - (c + c k + 2)y + + 1 \ / (c - v + l)(-c fe + v - 1) 

12 V (c + v + l)(c k + v + l) ) + 12 V (c + c k + 2)(y a -l) 

-2Li 2 ((c-v + W-J + v-l)) _ 2Li2 ( o-v + \) _ 2Li2 (* + V + + 1 



c + c fe + 2)(y+-l) 7 V c + y+ + V " Vc + f + 1 

2 log ( v + y A bg ^ (c + y a + l)(c fc -2/ Q + l) 



v-y a J \-y a {c + c k + 2) + cc k + c + c k + v 2 + 1 
I log f (c + y+ + l)( Cfc - y+ + l 



v -y+J \-y+( c + c fc + 2) + cc fc + c + c fc + v 2 + 1, 
+(log(l - j/ a ) - log(l - y+))(-2 log((c + v + l)(c fc + v + 1)) + 2 log(c + c fc + 2) 

+ log(l - //,,) + log(l - y + )) + 2Li 2 ( C ~ V + \ ) + 2Li 2 ( + // " r 1 



, c + y Q + 1 / \ c + i> + 1 

-(log(c + y a + l)- log(c + y+ + 1))(2 log(c + u + 1) - log(c + y a + 1) - log(c + y+ + 1)) 
+ (log(c fc - y Q + 1) - log(c fc - y+ + l))(21og(c fe + v + 1) - log(c fc - y a + 1) - log(c fc - y+ + 1)) 

-2Li 2 ( Cfc "" + M - 2Li 2 f Cfc -^ + 1 "l + 2Li 2 ( Cfc -" + M + 2Li 2 ^ " y + + 1 



(A.21) 



with 



y a = V (c fc (l - a) + (y 7 ^ - 1) a + l) 2 - c 2 + a(c fc - ^ + 1) . (A.22) 



For a massless emitter and a massive spectator the correction term is given by 

A/ Cik (0,^, a) = log(l - y a ) log ^- -^—A- — j log(y a + 1) log (- (fi 2 - l) (y a + 1)) 

+ 2Li 2 (—^—) ~ log(2 - 2^ fc ) log f — ) + log f log (— - 1 

- 2Li 2 (y*±l) , (A.23) 
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with 

ya = a(l - fj, k ) + V(a - l) ((//fc - l) 2 « - gfc + l) 2 ) (A 24) 

Mfe + 1 

Finally for the case of a massive emitter and a massless spectator we find 
A/ eik (/x i ,0, a) = -Li 2 - j^J - Li 2 (-a/i 2 + m| + a) + Li 2 ^1 - 

7T 2 

- log ((/i 2 - 1) (a - 1)) log (M 2 (-a) + m| + a) + 21og(^) log(a) + —. (A.25) 

In the same way as for the eikonal integral we take the a-dependence of the collinear 
integral into account by writing 

Z$' fc (Mj, e, ") = MfcJ e) + tefflkfaj, Mfc! a). (A.26) 

In the general case where both emitter and spectator are massive the correction term reads 
as 

2(q - l) Cj y + + ( Cj + 2y + )(2ay + + Cj ) (log + 41og(a)) 

AJ^.^a) = 2(c J+ 2y + )(2ay + + Cj ) ' 

(A.27) 

The limits \ij —> and Hk — > can be performed easily and one finds 

(a - l)^ 2 - a 



]_ (a-1) [nj-ljuj 
AI°f k ( H ,0; a) = -- t - _ \ 2_ lo S (~ a ^j +a + A*?) + 41og(a) V , 



and 



(A.28) 



Ai^i(0,// fc ;a) = --log(a) (A.29) 



respectively. 

For these calculations we have used the HypExp package for Mathematica [71,72]. 
As mentioned above the a-terms can as well be calculated using the formalism of mass 
regularization. To check our result we performed a second calculation using mass regular- 
ization and making use of the formulae giving in the appendix of [30] we found agreement 
between the two approaches. 
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